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ABSTRACT 

Aims. The stochastic acceleration of subrelativistic electrons from a background plasma is studied in order to find a possible 
explanation of the hard X-ray emission detected from the Coma cluster. 

Methods. We calculate the necessary energy supply as a function of the plasma temperature and of the electron energy, and 
we show that, for the same value of the hard X-ray flux, the energy supply changes gradually from its high value for the case 
when emitting particle are non-thermal to lower values when the electrons are thermal. The kinetic equations we use include 
terms describing particle thermalization as well as momentum diffusion due to the Fermi II acceleration. 

Results. We show that the temporal evolution of the particle distribution function has, at its final stationary stage, a rather 
specific form. This distribution function cannot be described by simple exponential or power-law expressions. A broad transfer 
region is formed by Coulomb collisions at energies between the Maxwellian and power-law parts of the distribution functions. In 
this region the radiative lifetime of a single quasi-thermal electron differs greatly from the lifetime of the distribution function 
as a whole. For a plasma temperature of 8 keV, the particles emitting bremsstrahlung at 20 — 80 keV lie in this quasi-thermal 
regime. We show that the energy supply required by quasi-thermal electrons to produce the observed hard X-ray flux from 
Coma is one or two orders of magnitude smaller than the value derived from the assumption of a nonthermal origin of the 
emitting particles. This result may solve the problem of rapid cluster overheating by nonthermal electrons raised by Petrosian 
(2001): while Petrosian's estimates are correct for nonthermal particles, they are inapplicable in the quasi-thermal range. We 
finally analyze the change in Coma's Sunyaev-Zeldovich effect caused by the implied distortions of the Maxwellian spectrum of 
electrons, and we show that evidence for the acceleration of subrelativistic electrons can, in principle, be derived from detailed 
spectral measurements. 
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1. Introduction 

An excess of hard X-ray (hereafter HXR) emission above 
the thermal spectrum has been found in the energy 
ran ge 20 — 80 keV from the Coma cluster of galax- 
ies JFusco-Femiano et aflll999llRephaeli fe Grubeill2002l 
iFusco-Femiano et all l2004|) . The v alidity of this excess 
is, ho wever, still unclear because iRossetti fe Molendil 
l)2004|) re-analyzed the Coma data and found no evidence 
for such h ard X-ray excess. Further r e-analysis of the 
same data l|Fusco-Femiano et all l2004|) seem to confirm 
the presence of the HXR excess at the level observed 
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bv IFusco-Femiano et all (|l999|) and iRephaeli fc Gruben 
(|2002h . The recent INTEGRAL observation of Coma 
(Renaud et al. 2006) analyzed the morphology of the clus- 
ter in the range 18-30 keV and found that the hard X-ray 
emission comes from an extended source with a radius 
~ 30 arcmin. The spatial distribution is similar to the 
thermal one as obtained with XMM (in the range 0.3-2 
keV). The INTEGRAL data indicate that the upper limit 
in 30-50 keV range is a factor w 1.5 above the mean RXTE 
spectrum and the non-thermal mechanisms are expected 
to contribute ~ 50% of the total flux in this region. 

The hard X-ray excess has been interpreted as 
being due cither to inverse Compton (IC) scatter- 
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ing of rel ic CMB photons by relativistic electrons 
(see, e.g.. ISarazin fc Lieul 11993 : iFusco-Femiano et all 
Il999t iBrunetti et all l200l|) or to bremsstrahlung of 
nonthermal subrelativistic electrons (|Enfilin et all Il999t 
iKempner fc Sarazinl l200fj(l . Yet further interpretations 
have b een proposed — X-ray emissio n by secondary elec- 
trons i Blasi fc Colafrancescd Il999|. or bremss trahlung 
emission by subrelativistic protons ijDogielL l200ll) . 

Each of these models has serious problems. In the 
framework of the IC model the magnetic field strength 
can be estimated from the ratio of X-ray to radio fluxes 
because both radiations are produced by the same rel- 
ativistic electrons. A weak (uniform) value of the mag- 
netic field of order 0.1-0.2 fiG is derived from the IC 
model. On the other hand, estimates of the magnetic 
field strength in the intracluster medium determined from 
Faraday rotation yield much higher value o f the order of 
~ 5—10 uG in the cluster ce nter (see, e.g.. lClarke et all 
l200lUOovoni fc Ferettl l2004h . However, modifications of 
the IC model such as the c omplex electron spectrum model 
ijSchlickeiser et al.1 Il987|) or the model with anisotropic 
pitch angle distribution of emitting particles I PetrosianL 
l200l|) may modify the estimates of the IC model quoted 
above. 

Models where the X-rays are generated by nonthermal 
electron and proton brcmsstrahlung are associated with 
an unacceptably large energy o utput of emitting particles 
I DogielL EoOll IPetrosianl Effll. 

However, the bremsstrahlung model has not been com- 
pletely explored because emission can also be produced by 
the quasi-thermal electron component which arises natu- 
rally when particles are accelerated from the background 
thermal plasma. In this case a part of the spectrum is 
formed under the influence of both Coul omb collisions an d 
a run-away flux of accelerated particle s llCurevichlll960ll . 
This class of mode ls was developed bv iDogiell (|2000l) and 
iLiang et alJ (|2002() , who assumed that the hard X-ray flux 
from a galaxy cluster is generated in regions of electron 
in-situ acceleration from the thermal pool. In such mod- 
els the electron distribution function develops an extended 
transition populated by quasi-thermal electrons. In this re- 
gion, which lies between the thermal and nonthermal parts 
of the spectrum, the distribution function differs strongly 
from the Maxwellian form because it is not an equilib- 
rium distribution (and in this sense it is not thermal), but 
is formed entirely by the Coulomb collisions (and therefore 
we cannot define it as nonthermal). 

The bremsstrahlung output from the emitting parti- 
cles is proportional to the lifetime of their distribution 
function. For electrons with totally nonthermal energies, 
this lifetime is about the timescale on which a single non- 
thermal electron suffers ionization losses. In this regime 
Pctrosian's (2001) estimates of the energy output required 
to reproduce the hard X-ray excess of Coma is completely 
correct. However, for a given energy in the regime where 
the spectrum is formed by Coulomb collisions, the life- 
times of the distribution function and that of a single elec- 
tron may differ dramatically. In fact, the lifetime of, e.g., 



an equilibrium spectrum is much larger than the char- 
acteristic lifetime of a single particle of this spectrum. 
Therefore, if the hard X-ray Coma flux is emitted by elec- 
trons from the quasi-thermal transfer region, a specific and 
more detailed analysis is necessary for such a situation. 

As it was shown by IDogiell l)2000(l and ILiang et alJ 

i|2002|) . if the plasma temperature is of order of several 
kcV, then hard X-rays in the energy range 20 — 80 keV 
are produced by this quasi-thermal component of the elec- 
tron flux. For a fixed radiated luminosity, such as that of 
the hard X-ray flux from Coma, the rate of energy input 
into the electrons can be lower than that required in non- 
thermal bremsstrahlung models. We shall discuss below 
this energetic problem in more detail. We will show, in 
addition, that this acceleration process can be tested by 
looking at the spectral changes in the associated Sunyaev- 
Zcldovich (SZ) effect from the intracluster medium. 

The problem here described will be analyzed in the 
present paper under the following assumptions: 



1. The hard X-ray flux is produced by bremsstrahlung of 
subrelativistic electrons; 

2. Electrons are in-situ accelerated from background 
plasma by a stochastic, Fermi-II type acceleration 
mechanism; 

3. The characteristic time of acceleration is much larger 
than the characteristic time of Coulomb collisions of 
thermal particles, i.e. a small part of background par- 
ticles is accelerated. Therefore, the parameters describ- 
ing the background plasma change very slowly; 

4. This allows to neglect the nonlinear terms in kinetic 
equations; 

5. The particle acceleration is investigated in subrelativis- 
tic energy range. We do not consider, therefore, accel- 
eration of fully relativistic electrons. 



The first part of the paper is focused on the energy 
problem. We will show that the excess of hard X-ray emis- 
sion above the thermal spectrum of Coma can be produced 
by electron bremsstrahlung and that the necessary energy 
output of electrons can be smaller than in the nonthermal 
bremsstrahlung model. Our analysis is based on simplified 
analytical solutions of the kinetic equations. Nonetheless, 
beyond the simplified description of the problem, it shows 
that the crucial energy problem set by the HXR excess 
can be definitely alleviated. 

In the second part of the paper we derive the charac- 
teristics of in-situ acceleration from the observed spectrum 
of X-ray emission and we calculate the distortions of the 
equilibrium Maxwellian spectrum due to this acceleration. 
We finally explore whether a signal of acceleration in the 
Coma halo can be indicated by the associated Sunyaev- 
Zcldovich effect. In order to compare our results directly 
with those of Petrosian (2001), we use here the same cos- 
mological model with Hq — 60 km s _1 Mpc" 1 . 
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2. Kinetic Equation for Electrons 



The dimensionless rate of ionization loss has the form 



Shock acceleration is usually considered as a candidate 
for particle productio n in the intracluster medium (see, 
e.g., IKuo et all 12003^ . and for some clusters, the avail- 
able observations suggest that there are shocks (see, 
e.g., iMarkevitch et al. L 120051 iFabian et all. l2006|) and 
elect rons which mi g ht be accelerated by shoc k s (see , 
e.g. iRrunetti et all 120011 iMiniati et all 120011 1200,^ . 
However, extended radio and (perhaps) X-ray emission 
cannot be associated with strong shocks since accelerated 
electrons are unable to travel large distances from their 
sources without loosing much of their energy. 

Plasma turbulence is hence considered as a vi- 
able model for particle (re) acceleration in cluster halos. 
Numerical calculations sh ow that strong turbulenc e can be 
excited in halos (sec, e.g., iRicker fe SarazmLl200l|) . If this 
turbulence generates plasma waves, then a slow stochastic 
process can accelerate particles through their resonant in- 
teractions with the waves. This process can be described 
as momentum diffusion. We note that particles can also be 
accelerated directly by hydrodynamic turbulent or quasi- 
periodic flows in a manner simi lar to stochastic ac celera- 
tion by plasma waves (see, e.g.. lWebb et all 12003^1 . 

The evolution of the distribution function of particles 
which are scattered by electromagnetic fluctuations is de- 
scribed by the Fokkcr-Planck equation which can be trans- 
formed to the diffusion type equation by integration over 
particle pitch-angles, if scattering is very effective and the 
distribution function is quasi-isotropic. For the mechanism 
of the in-situ acceleration from background plasma, the 
equation can be written in the form 



df 
dt 



i__a_ 

p 2 dp^ 



{D c (p)+D(p)} 



dj_ 

dp 



0,(1) 



where (dp/dt)c and D c (j>) describe particle convection 
and diffusion in the momentum space due to Coulomb 
collisions, and D(p) is the diffusion coefficient due to the 
stochastic acceleration. Eq. Q can be written, in a dimen- 
sionless form, as 



df 

dt 



dp 



B(p).f 



0. 



(2) 



where p — p/\/mkT is the dimensionless momentum, t = 
tu is the dimensionless time and D p (p) = D p {p) / (vmkT) 
is the diffusion coefficient. The frequency v is 



1"Knc rim 



(3) 



V mkT 

where r e = e 2 /(mc 2 ) is the classical electron radius. Here 



B(p) = p 



and 



~ 2 / d P 



dt 



A{p) = B(p) 



kT 



P 2 D p (p) 



(4) 



(5) 



V 2 



d£\ = 1 
dt) i P 

Emc 2 ^ 2 



kT 



(6) 



In 



h 2 io 2 Y 



0.43 



where uj p is plasma frequency and E(p) is the particle 
total energ y The quasi-stea dy state solution of Eq. (0), 
obtained bv lr T urevichllll9fi0h . reads 



/ = \j-n(t)exp 



B(v) 
A(v) 



dv G(p) 



(7) 



where 



G(p) 



J[dv/A(v)]e X p [J[B(t)/A(t)]dt 

OO / V 

J[dv/A(v)]e X p(j[B(t)/A(t)}dt 
o \o 



(8) 



and n(t) describes slow variations of the background 
plasma density, consistently with the assumption of slow 
acceleration. 

Detailed information on the conditions necessary to de- 
rive a reliable momentum diffusion coefficient D(p) are not 
well determined yet. First of all, its value is determined by 
a spectrum of electromagnetic fluctuations W(k) (where 
k is the wave-number of fluctuations) which is basically 
unknown, t hough new th eoretical treatments of plasma 
turbulen ce llVermal 1 19991) or X- ray observations of clus- 
ters (see ISchuecker et al. 1 12004|) indicate a Kolmogorov- 
Oboukhoh type spectrum of turbulence in the large scale 
range, between 20 kpc and 2.8 Mpc. Secondly, we do not 
know in details the ratio between the energy density of 
the intracluster plasma, Wth — 1 eV cm -3 , and that of 
the magnetic field, Uh- This ratio is, moreover, totally 
unknown in the regions of particle acceleration. For the 
magnetic field strength observed in clusters (whose esti- 
mates ranges from ~ 0.1 to ~ 10 /j,G, (see e.g. Carilli 
& Taylor (2002) for a review), the ratio (3 = W th /U H 
is in the range ~ 0.4 (low (3 plasma) up to ~ 1000 
(high (3 plasma) . Additional theoretical uncertainties come 
from the lack of knowledge of the mechanism through 
which turbulence is formed in regions of particle accel- 
eration, whether it is developed by cascade processes 
(as it is in the new mode l of turbulence developed by 
iGoldreich fc Sridharl lll997l)). bv intermed iate turbulence 
model (see, also ICho and Lazarianl . l2004h or if it is due 
to the flux instability when a flux of particles escap- 
ing from acceleration regions excites there MHD fluctu- 
ations due to reso nant interaction, a s it may occur in the 
Galactic halo (see iDogiel et allll994|) or near shock fronts 
l(Ptuskin fc Zirakashvili LI200!tI) . Given all these uncertain- 
ties, we are unable to choose reliable parameters of the 
kinetic equations, and we have to resort to a quite gen- 
eral description whose overall features can be, nonetheless, 
tested against the available data. 
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In the simplest case of charge particle scattering, the 
momentum diffusion coefficient has the form 



D p (p) = D p 



(9) 



for nonrelativistic particles (see iToptvgirl Il985|) . 

For the case of low-/? plasma the momentum dif- 
fusion coefficient was der i ved b y iMiller fc Steinackerl 
<ll992l): ISteinacker fc Millerl lll992l) : IMiller et all <ll99filt : 
I S chlickers^r*^T^iUerl ~ i|l99q) for resonant particle-wave 
acceleration in solar corona and in the interstellar 



medium of the Galaxy. It s analytical form taken from 
ISchlickeiser fc Millerl ill 99^ is 



D P (P) 



7r(j-l)ci(a) fm 
4 \U H 

U{r L k min y — 5- In — , 



(10) 



where tl and SI are the Larmor radius and cyclotron fre- 
quency of protons, k m i n is the minimum wave number of 
the MHD spectrum, W t is the total energy of the magnetic 
fluctuations and Ur is the energy density of the large scale 
magnetic field. Here the spectrum of magnetic fluctuations 
is supposed to be described by a power-law, W(k) oc k~i , 
where ci(j) is a constant depending on the spectral index 
J- 

The momentum diffusion coefficient in a high-/? plasma 
is determined by the spectrum of magnetic turbulence ex- 
cited by stochastic plasma motion. In the approximation 
of strong turbulence, corresponding to high-/? plasma, its 
value depends on whether the accelerated particles are 
magnetized or unmagnctized in a random magnetic field. 
The coefficients of kinetic equations are determined by 
pair-correlations of random velocity and random magnetic 
fluctuations. Note that there is no resonant interaction 
in this case. Below we present general equa tions for the 
mome ntum diffusion coefficient derived by iDogiel et alJ 
l)l987j) . For magnetized particles the coefficient of momen- 
tum diffusion is given by 



From an estimate of the acceleration time scale, we can 
derive the energy density of the resonant magnetic fluc- 
tuations in the case of low-/? plasma and the correlation 
length of turbulence in the case of high-/? plasma. 
At the present stage of our knowledge, we cannot prove 
or disprove the validity of any form of the diffusio n coef- 
ficient; however, we notice that lClarke et all 1 2f)0ll) found 
from ROSAT and radio observations that the total mag- 
netic energy content in clusters is comparable to the total 
thermal energy content in the same cluster volume. 

To circumvent this problem, we derive here the gen- 
eral characteristics of the acceleration mechanism - i.e., 
the characteristic time of acceleration (the dimensionless 
parameter a, see below) necessary to produce the X-ray 
excess above the thermal distribution - from the observed 
flux of hard X-rays from Coma, a procedure which is in- 
dependent of the details of the acceleration process. 

In the following, the spectrum of turbulence is as- 
sumed to be a power-law, W(k) oc fc"- 7 , with j ~ 
2, i.e. a spectrum th at is close to that derived by 
ISchuecker et al. 1 1)2004]) . For this fluctuation spectrum the 
momentum diffusion coefficient D{p) is also a power-law 
function of p in a confined range of momenta as fol- 
lows fr om IMiller fc Steinackerl l| 19921) ; ISteinacker fc Millerl 
(|l992lh see.Fig. 5). 
The energy of accelerated particles 



W e OC / p 4 f(p)d P: 



(14) 



diverges for any reasonable power-law form of the diffu- 
sion coefficient as follows from Eq.(J2J). In order to avoid 
this divergence we should make the natural assumption 
that there is a cut-off in the fluctuation spectrum for a 
wave number fc that gives the maximum momentum p of 
accelerated particles. We discuss the expected solution us- 
ing a simple form of the momentum diffusion coefficient 



£>_ = 6 / < VV' >< hti >< VV'hh' > dt' 



\ 0^ 



p <p 
p>p 



(15) 



where V is the turbulent velocity of the plasma, h = H/ | 
H | is the random direction of the magnetic field line. The 
values of < VV >, < hh' >, and < W'hh' > arc the 
pair-correlations of turbulent velocity, direction of mag- 
netic field and its derivative, respectively. For unmagne- 
tized particles one obtains 



D„ 



2e2 
3c 2 



< V >t.vt< H >t.vt dr . 



(12) 



where 

< V 2 > T ,vT= I | V(k, w) | 2 cxp (i(w - kv)r) dud 3 k .(13) 



The boundary conditions at p = p are discussed in the 
Appendix. 

We set the cut-off momentum p at an energy E = 
p 2 /2m > 80 keV, above the range of observations, so that 
it does not affect the calculated bremsstrahhmg spectrum 
. In our calculations we take the minimum acceptable 
value of the cut-off energy to be Ek = 110 keV. In the 
case of resonant interaction this energy corresponds to the 
minimum wave-number fc m .;„ ~ 4.3 • 10~ 9 cm -1 in the 



1 As follows fro m observations of t he X-ray flux from the 
Galactic disk (see lLebrun et ail 12004) , similar processes may 
accelerate particles up to ^ 100 keV in the interstellar medium, 
although one cannot exclude the possibility that electrons are 
accelerated to far higher energies. 
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spectrum of magnetic fluctuations W(k). Therefore, the 
energy density of resonant waves is 



Wt 



W(k)dk. 



(16) 



In the nonrelativistic energy range, eq.© is simplified sig- 
nificantly and takes the form 



di 



1_ d_ 

p 2 dp 



1 



b{p)p< 



of 



dp 



'- + f 



o. 



(17) 



For the qualitative analysis of the problem we consider the 
dimcnsionless diffusion coefficient in the power-law form 
£>{p) = D{p) / {vimkT) = ap q that enable us to get simple 
analytical solutions of ea. (|17|l. Thus, for q = 1 the solution 
of cq. H17fl has a simple form ijGurevichl Il960|) 



f{p) oc exp 



Note that here 



Anne 



(fcT) 3 /2 



2 arctan( v / ap 
4V" 



In A, 



(18) 



(19) 



where In A is the Coulomb logarithm. 

The dimcnsionless parameter a f or Coma can b e de- 
rived from observational data. Thus, iDoeiell l)2000|) ana- 
lyzed electron acceleration in the ce ntral part of Coma 
with the density n ~ 10 -3 cm -3 while llTang et all l|2002fl 
derived parameters of acceleration for the average gas den- 
sity in the Coma halo, n ~ 10 -4 cm -3 . It is rather difficult 
to compare results of these investigations because different 
forms of momentum diffusion were used in these publica- 
tions. 

We notice that the solution (|18fl provides an illustra- 
tive oversimplification of the solution of Eqs.Q which is 
useful in order to get rough quantitative estimates. We 
present in the following numerical calculations of Eq. (|17fl 
which show the time variation of the distribution func- 
tion / under the influence of stochastic acceleration and 
Coulomb collisions. 



3. Time Variations of the Spectrum of 
Accelerated Particles 

We solve Eq. (|17fl numerically in order to understand the 
time evolution of the distribution function / for a diffu- 
sion coefficient D(p) = ap q with a w 0.001 and q = 1. As 
discussed earlier, D(p) must have a cutoff for large enough 
momentum. For simplicity, we choose the cutoff momen- 
tum to be larger than the maximum momentum we used 
in our numerical calculations. In fact, the position of the 
cutoff momentum does not significantly affect the evolu- 
tion of the distribution function at small momenta. 

We considered two initial cases: i) transformation 
of the Maxwellian distribution under the influence of 
stochastic acceleration. In this case we can estimate the 
characteristic time for the formation of the particle excess 



above the thermal distribution; ii) transformation of the 
resulting noncquilibrium spectrum under the influence of 
Coulomb collisions only. 

These calculations allow us to estimate the power of the 
electron source necessary to compensate collisional dissi- 
pation and to keep the particle excess above the thermal 
distribution at the level necessary for the production of 
the observed hard X-ray emission in Coma. 

In both cases, we calculate the evolution of the electron 
spectrum up to the dimensionless time i = 4000 (which is 



equal to 1.85 Gyr for Coma, where v\ f=a 7.2 x 10 



As for the boundary conditions, we use a free boundary 
condition at the high-momentum boundary of our calcu- 
lations, while at the low momentum boundary we used 
zero-flux boundary conditions. 

Figure 1 shows the temporal evolution of the distribu- 
tion function / formed by the combined effects of stochas- 
tic acceleration and thcrmalization due to the Coulomb 
collisions. The distribution function is normalized to the 
value at the low momentum boundary. The dashed curve 
shown in this figure is the original Maxwellian distribu- 
tion, and the dotted curve is the steady state solution 
(with a cutoff at large momentum). The five solid curves 
represent the distribution function at increasing times 
t = 800, 1600, 2400, 3200, and 4000. The figure shows that 
the distribution function rapidly approaches the steady 
state solution for the relatively low momenta of the quasi- 
thermal particle regime (j> 10). However, for nonther- 
mal rclativistic particles the stationary state is reached on 
timescales longer than a cluster lifetime. Therefore, in the 
framework of this model it is difficult to expect that sub- 
relativistic electrons emitting hard X-rays and relativis- 
tic electrons emitting radio emission are produced by a 
single mechanism of particle acceleration. Therefore, this 
analysis cannot be extended to relativistic energies. Our 
assumption that the maximum energy of the accelerated 
electrons is 110 kcV is hence justifiable on the basis of the 
previo us results. Moreover, as shown by IWolfe fc MeiiH 
( 2006) based on a covariant kinetic theory of electron plas- 
mas, a power-law tail obtainable from direct or stochas- 
tic acceleration of relativistic particles cannot survive for 
times longer than ~ 20 Myrs because the equilibration 
time scale for relativistic electrons is quite short for the 
case of clusters and hence a thermal distribution is soon 
established. 

Now let us assume that the acceleration mechanism 
is interrupted but the initial distribution function is de- 
scribed by the noncquilibrium form of ea. l|18|) . In this case, 
the dimcnsionless equation for the distribution function / 



di 



1__5 

p 2 dp 



ld_l + f 
p dp 



(20) 



Figure 2 shows the evolution of the electron distribu- 
tion function, as described by Eq. J3UJ), if acceleration is 
switched off and the distribution function is allowed to 
evolve from its quasi-steady state form under the influ- 
ence of collisional dissipation. The collisional regime of 
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Fig. 1. Temporal evolution of the distribution function up 
to t = 4000 (1.85 Gyr) after acceleration is switched on. 
The initial distribution is Maxwellian. 



Fig. 2. Temporal evolution of the distribution function 
up to t = 4000 (1.85 Gyr) after acceleration is switched 
off. The solution (|18|l of Eq. (J2J) was taken as an initial 
distribution function. 



quasi-thermal particles lies in the range 5.5 < p < 30 for 
a = 0.001 and q = 1. 

We show in Fig. 2 that the dissipation time scale is larger 
than the single electron ionization loss time scale near the 
thermal particle region, because the distribution function 
there is of almost equilibrium form, and it increases in 
the region of nonthermal momenta, where the lifetime of 
particles increases as Ti oc p 3 . 

4. The Energetics of Quasi-Thermal Electrons 
and the Coma hard X-ray Flux 

We now return to the problem of the energetics of the 
emitting electrons. First, we recall Petrosian ' s crit icism 
of the bremsstrahlung model (see IPetrosianL l200l|) . He 
estimated the yield of bremsstrahlung photons as 7 ~ 
(dE/dt)br/(dE/dt)i - 3 x 10~ 6 in the energy range 20-80 
keV, where (dE/dt)i/ (dE/dt)b r is the ratio of ionization to 
bremsstrahlung losses. Then for the hard X-ray flux from 
Coma in this energy range, F x ~ 4x 10 43 erg s _1 , a large 
amount of energy, F e ~ F x /Y ~ 10 49 erg s _1 is trans- 
ferred from the accelerated electrons to the background 
plasma by ionization losses. As a result, the intracluster 
plasma temperature should rise to a temperature > 10 s K 
on a quite short time scale ~ 3 x 10 7 yrs. We stress here 
that these conclusions were obtained under the assump- 
tion that the lifetime of a single electron equals the lifetime 
of the particle distribution function. These estimates are 
correct only in the case that the electrons are nonther- 
mal, and therefore collisionless. However, the previous en- 
ergetic arguments cannot be used in energy ranges where 
the spectrum is formed by Coulomb collisions because, as 
we have shown in Scction[j2 the lifetime of particles differs 



strongly from the lifetime of the distribution function (see 
FigQJ. For instance, let us consider the lifetime of thermal 
electrons at an energy of about 8 keV. Their individual 
lifetime is about ~ 10 6 yr, but the lifetime of the distri- 
bution function at these energies is much longer (almost 
infinite) because the distribution function for these ener- 
gies is almost in equilibrium. Therefore, the energy supply 
necessary for the bremsstrahlung radiation can only be es- 
timated from the corresponding kinetic equations. It fol- 
lows that estimates of the energetics based on the lifetime 
of single electrons are not appropriate here, and lead to 
wrong conclusions. Figure 3 shows, in fact, the variation 
of the lifetime of the distribution function in eq. Ijl8|l . as 
derived from Eq. (|20|l . From this figure we see that the 
lifetimes of the distribution function and of the particles 
are equal to each other only for high (nonthermal) electron 
energies. 

In order to estimate the energy supply necessary 
to support the non-equilibrium distribution fca ll8l) we 
use the following kinetic equations in which we include 
bremsstrahlung losses 



dW e _ 



= $ 



dt 

= AttVE, 



(21) 



15/ 

p dp 



l+p 2 



br 



f 



Here the total number of particles with momentum > p 
in a volume V is F e {p) = V J p f(p)4Trp 2 dp, and the total 
electron energy in this volume is W e ~ E e F e , where the 
energy E e = kTp 2 /2. Then from Eq. (|2Tj|l we obtain an 
expression for the rate of change of the energy content of 
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Fig. 3. The lifetime of the electron spectrum l|18|) derived 
from Eq. (|17J) is shown by the solid curve. The lifetime of 
a single electron determined by ionization loss is shown by 
the dashed-dotted curve. Here the dimension frequency of 
Coulomb collisions of thermal particles v = 7.2- 10~ 14 s" 1 
for the Coma parameters T = 8.2 keV and n = 1.23 ■ 10 
cm -3 . 



the electrons. The number of electrons emitting 50 keV 
radiation is fixed in order to satisfy the observations. 

In the case of nonthermal electrons whose spectrum is, 
e.g., a power-law, the last term exceeds the others on the 
right-hand side of Eq. I|21[) . and then 



dW e 
dt 



n ! It 



$ ~4TrVE e f e -, 

Ti 



(22) 



which is the result discussed by Petrosian (2001), where 
the electron energy loss rate is determined by the ioniza- 
tion loss of nonthermal particles. Here r,; is the character- 
istic time scale of the ionization loss at the energy E e . 

If the particles are thermal, and their spectrum is de- 
scribed by a Maxwellian, then the first and the third terms 
on the right-hand side of eq.(20) cancel out so that 



dt 



th 



dj_ 

br d P 



We 
Tbr 



$0 



(dE/dt) t 
{dE/dt) br 



(23) 



Here, the time-scale Tb r is the characteristic time for 
bremsstrahlung loss, and Tb r 3> n. Hence, (dW e /dt) nt h S> 
(dW e /dt)th for the same flux of bremsstrahlung radiation 
produced by these electrons. 

We can estimate the rate of energy supply to the 50- 
keV electrons that is required to generate the observed flux 
of ~50-keV bremsstrahlung X-ray emission from Coma as 
a function of different values of the plasma temperature 



T . For simplicity, we consider the electron spectrum to be 
of the form given by eq. lfT%|) . When the temperature T is 
low, the 50-keV electrons are in the nonthermal particle 
regime, while for high plasma temperature these electrons 
are thermal. From Eqs. (|18f) and (|21|l we obtain the vari- 
ation of the associated rate of electron heating as a func- 
tion of the background temperature T . The required level 
of heating, normalized to the rate of heating required for 
non-thermal particles, <&o = (dW e /dt) n th, is shown as a 
function of the temperature in Fig. where the parame- 
ters defining the acceleration of particles are those inferred 
from the X-ray analysis of the Coma cluster. We find that 




0.005 



T(keV) 

Fig. 4. The rate at which energy must be supplied to the 
50-kcV electrons generating a fixed bremsstrahlung X-ray 
flux as a function of temperature of background gas, nor- 
malized to $o, the rate at which energy must be sup- 
plied if the emitting electrons are non-thermal. For Coma 



$o ~ 10 49 erg s" 1 l(PetrosianLl200l 



when the temperature T is low, the 50-keV electrons are 
nonthermal and - as expected - the rate at which they 
must be heated is independent of the temperature and 
is almost equal to $o- However, when the plasma tem- 
perature T increases, and the 50-keV electrons shift into 
the quasi-thermal regime (where the spectrum is formed 
by the Coulomb collisions) , the required energy supply to 
maintain the observed hard X-ray flux from Coma de- 
creases rapidly as the background plasma temperature 
increases. For a temperature T ~ 8 keV, the emitting 
electrons in Coma are quasi-thermal and the energy sup- 
ply they require is one or two orders of magnitude below 
Petrosian's (2001) result, even by using a simple qualita- 
tive estimate. Accurate quantitative calculations may give 
even higher variations. This result may indeed solve the 
energetic problem raised by Petrosian (2001): in fact, the 
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quasi-thermal electrons need much less energy in order to 
produce the observed HXR bremsstrahlung radiation with 
respect to the case of nonthermal particles. 

As it is clear from our analysis, our model describes 
the processe of in-situ acceleration reasonably well for 
relatively long time because of the presence of weak ac- 
celeration mechanisms. Therefore, the time variations of 
both the plasma density and of its temperature are very 
slow. Attempts to investigate a nonlinear phenomenolog- 
ical model o f particle acceleration in Coma were made by 
iBlasi I l|200d|) . Such an analysis showed that, due to nonlin- 
ear processes, the temperature of plasma increases slowly. 
However, there arc still questions on whether this model 
can describe correctly the process of ther malization in the 
cluste r atmospheres (see discussion by IWolfe fc Melia L 
l200fij) . The role of non- linear effects certainly needs fur- 
ther analysis which goes beyond the scope of this paper 
and we will address this issue elsewhere. 




5. Bremsstrahlung emission of quasi-thermal 
electrons 

We derive here the spectrum of the emitting electrons in 
Coma from the X-ray data using the more accurate solu- 
tion (see eqCJ) of the kinetic equation eq.(0). We calculate 
the flux of hard X-ray emission using the equation 

oo 

F "^fJ n 'k N{E)iE - (24) 

with cross-section 

do x 16 e 2 rime 2 
dE^ ~ Yhc 1 EE X 

xb p^ ), (25, 

The calculated bremsstrahlung spectrum from Coma 
is shown together with the Beppo-SAX data in Fig. 
We consider a the volume of particle acceleration to be 
V ~ 7.74 x 10 74 cm 3 which is smaller than the whole vol- 
ume of the Coma halo inferred from radio and soft X-ray 
data. In the general case it is convenient using the simplest 
form of the diffusion coefficient in cq.@. The necessary 
acceleration time of ~ keV particles, r acc = p 2 /D(p), is 
of the order ~ 10 17 s, independently of the acceleration 
m echanism. Thi s value is slightly lower than that derived 
bv lDogiel (|2000Tl for the central region of Coma where the 
intracluster gas has the highest density. 

For the case of low-/3 plasma we can go a step 
further and estimate the necessary density of resonant 
wave by using the diff usion coefficients for elec trons 
and protons taken from | Steinacker fc Milled l|l992l) and 
ISchlickeiser fc Milled l|l99a) . We notice that the momen- 
tum diffusion coefficient derived from these equations cor- 
responds to the general form © in the nonrelativistic en- 
ergy range. 



Fig. 5. Bremsstrahlung X-ray flux from the Coma cluster 
produced by subrelativistic electrons with a cut-off, E = 
110 keV, T = 8.25 k eV. The flux measured with Beppo- 
SAX was taken from iFusco-Femiano et alJ l|l999j) . 
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Fig. 6. The dimcnsionlcss momentum diffusion coefficient 
D(p) = D(jp) / (ymkT) is shown as a function of the di- 
mensionless momentum p = p/ V mkT. 



Fig. 6 shows the momentum diffusion coefficient for elec- 
trons for the low-/3 case. Only a low energy density in res- 
onant plasma waves, (Wt/Un) — 3.8 • 10~ 13 , is required 
to fit the Coma HXR data. Here Wt is the energy density 
of plasma waves and Ur = Hq/8it is the energy density 
of the large scale magnetic field. 
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The corresponding electron distribution function is 
compared with a thermal spectrum at T = 8.25 keV 
in Fig. 7. The change in the distribution function is 
achieved by ~ 10% of the background electrons being 
quasi-thermal. It follows that the energy density of ther- 
mal electrons is W thermal ~ 1.6 eV cm~ 3 and that of the 
quasi-thermal electrons is Wnonthermal ~ 0.56 cV cm -3 , 
thus yielding the ratio W nonthermal I W thermal ~ 0.35 for 
Coma. 



0.00001 
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Fig. 7. Spectrum of subrelativistic electrons with cut-off, 
p = 110 keV, and T = 8.25 keV derived from the X- 
ray flux from the Coma cluster (solid curve) . The thermal 
spectrum is shown by the dashed-dotted curve. 



~ 8 x 10 48 erg s 1 , which matches the requi red rate if pro - 
tons produce the observed hard X-ray flux ijDogiell . l20mh . 
However, this process involves subrelativistic protons and 
so it faces the same problem of energetics and plasma heat- 
ing as we discussed in Section 0] 

The theoretical approach we have taken in the present 
paper allows us to estimate the density of electrons and 
protons accelerated from the background pool using the 
kinetic parameters derived from the spectrum of plasma 
waves, and hence to understand whether protons are im- 
portant. The momentum diffusion coefficient for protons 
was taken in the form of Ea. (flop . In Fig. 8 the spec- 




E(keV) 

Fig. 8. The spectrum of subrelativistic electrons (dashed- 
dotted curve) and protons (solid curve) in the Coma clus- 
ter halo are shown as a function of the particle energy. 



6. The spectrum of accelerated protons 

Protons are generally disregarded in the calculation of 
bremsstrahhmg emission because, having low rates of en- 
ergy loss, they leave the emitting region before they lose 
a significant fraction of their energy. On the other hand, 
if protons escape from a radiating region relatively slowly, 
they can produce a bremsstrahhmg flux whose value is 
comparable with that of the electrons. An analysis of the 
brcmsstrahlung radiation emitted by subrelativis tic pro- 
tons (inverse bremsstrahlung) in the Coma cluster ijDogielL 
l200lh showed that protons are able, in fact, to generate the 
observed hard X-ray emission. These X-rays are accompa- 
nied by the excitation of background nuclei which could 
be detectable through the develo pment of prom inent car- 
bon and oxygen gamma-ray lines. iDoeiell l|200ll estimated 
the ex pected flux of these lines. Recently. Ilvudin et alJ 
l)2004l) have found tracers of this gamma-ray line emis- 
sion towards the Coma and Virgo clusters at the ex- 
pected level. The e nergy deposit i on by subrelativistic pro- 
tons estimated by Ilvudin et all l|2004|) is of the order of 



trum of protons is shown for the parameters derived 
in Section [SJ It can be seen that stochastic accelera- 
tion is ineffective for subrelativistic protons (note that a 
similar conclusion fo r the Galaxy was also obtained by 
ISchlickeiser fc Milled (Il998h 1. In order to produce a flux 
of proton bremsstrahlung comparable with that of elec- 
trons, for example at a photon energy 10 keV, the density 
of 20-MeV protons should be about the same as the den- 
sity of 10-keV electrons. As shown Fig. 8, this condition is 
not fulfilled, and so we conclude that the hard-X-ray emis- 
sion from the Coma cluster cannot be ascribed to proton 
bremsstrahlung. 



7. The Sunyaev-Zeldovich effect produced by 
in-situ accelerated particles 

In the context of the present study it is important to 
find observational resources which can provide indepen- 
dent evidence for the subrelativistic electron population 
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that we consider to be responsible for the hard X-ray 
emission from clusters of galaxies. One possibility is to 
use a detailed analysis of the inverse Compton scattering 
of CMB photons off the population of sub-relativistic elec- 
tron s, the SZ effect described b y Zeldovich and Sunyaev 
( see IZeldovich fc Sunvaevl . Il969h . The amplitude and the 
spectrum of this effect depend on the distribution function 
of the electrons in the intracluster medium. As a result of 
the inverse Compton scattering, the spectrum of the CMB 
radiation is shifted to higher frequencies when observed 
along the line of sight t hrough the intrac luster medium 
(see for a general review lBirkinshawlll999l) . The spectral 
distortion of the microwave background can be calculated 
for any electron distribution function, and not only for 
the standard Maxwellian spectrum adopted in most dis- 
cussions of the Sunyaev- Zeldovich effect: hence, we can 
calculate the Sunyaev- Zeldovich effect as due to a nonther- 
mal component of the electron distribution in the supra- 
therm al (see, e.g.. lBlasi et all l200nUEnsslin and Kaizer L 
and relativistic (see IShimon and Rephaeli I l2002t 
Colafrancesco et al~l l2003|) energy ranges. To this aim, 
we use here the electron distribution function derived 
from the X-ray data to calculate the S unvae v- Z e l dovich 
effe ct using the forma lism presented in iRephaelJ l) 19951) 
and lBirkinshawl l|l999j) . The change in the radiation tem- 
perature AT(i>) at frequency v is given by 



AT(y) (e x - l) 2 AI(v) 



' rad 



(26) 



where x = hv/kT ra d is the dimensionless frequency, Al is 
the scattering-induced change in the specific intensity of 
the cosmic microwave background spectrum at frequency 
T ra d = 2.73 K is the temperature of the microwave 
background radiation, and the scale of the specific inten- 
sity is 



In = 



2(Wra<i) 

(he) 2 



(27) 



The CMB temperature variation due to the SZ effect is 



(e x - l) 2 

AT{v)=T rad T y > I dsP^s) 



(28) 



/ exp[— 3s)] 



\exp[xe( _s )] — 1 exp[x] — 1 

where r = otk J dln e is the optical depth along a line of 
sight of length I through electrons of density n e , with axh 
being is the Thomson cross-section, and 



Pe (f3)P(s,f3)d(3 



(29) 



is the photon redistribution function calculated in the 
limit of single scattering (appropriate here for low values 
of r), with 



A, 



M_l 



(30) 



The quantity p e ((3) is the normalized electron spectrum 
given as a function of the normalized velocity (3 = v/c, 



p e (/3)d(3 = 1 



(31) 



The function P{s,j3) describes the logarithmic frequency 
ratio cause d by a single electr on/photon scattering, and it 
is given in iBirkinshawl l|l999|) . The general description of 
the non-thermal SZ effect for multiple scattering, general 
va lues of r and multiple electr on distribution can be found 
in lColafrancesco et al. 1 l|2003() . For practical purposes, we 
transform the electron distribution function from its mo- 
mentum representation (Fig. 7) to its /3 representation 
(Fig. 9). The suprathermal excess of electrons which is 
evident in these figures is compensated by a reduction in 
the number of electrons at lower energies, but this reduc- 
tion is a small fractional change of the thermal electron 
number, and so it is neglected in Figs. 7 and 9. 




Fig. 9. The dimensionless distribution function for subrel- 
ativistic electrons with cut-off p = 110 keV and T = 8.25 
keV derived from the X-ray flux from the Coma cluster 
(solid curve). The pure thermal spectrum is shown by the 
dashed-dottcd curve. 



We can now calculate the Sunyaev- Zeldovich effect 
from the regions where the acceleration takes place. Here 
we should take into account that, as follows from the 
soft X-ray data, the volume of the acceleration region, 
Vacc = 7.7 x 10 74 cm 3 , is comparable with the Beppo-SAX 
estimate of the emitting volume V em . < 1.7 x 10 75 cm 3 . The 
total optical depth of t he thermal gas in the Co ma cluster 
is not precisely known. iBattistelli et all l|2003h estimated 
the total optical depth of Coma as r t h = (5.35 ± 0.67) x 
10~ 3 and iDe Petris et alJ (|2002() estimated the depth as 
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Fig. 10. The Sunyaev-Zeldovich effect produced by ther- 
mal electrons in the Coma cluster region, with temper- 
ature T = 8.25 keV and optical depth r = 5.3 x 10" 3 
(dashed curve). The solid curve shows the Sunyaev- 
Zeldovich effect due to thermal electrons from the Coma 
halo plus thermal and quasi-thermal electrons from re- 
gions of particle acceleration which occupy a small fraction 
of the Coma halo and have optical depth 9 x 10 -4 ). 



i~th = (4.1 ± 0.9) x 1CP 3 . Moreover, some contribution 
to the optical depth may be due to an extended halo of 
cooler gas, with a temper ature 0.6 — 1.3 keV and density 
of 10" 4 - 10~ 3 cm^ 3 Csee lNevalainen et allEoolj) . 

We calculate here the Sunyaev-Zeldovich effects from 
the thermal volume of the Coma halo in the absence of 
particle acceleration (the dashed curve in Fig. 10) and in 
its presence (the solid curve in Fig. 10). In these calcu- 
lations we assumed that the temperature and the optical 
depth of the thermal electron population arc T = 8.25 keV 
and Tth = 5.3 x 10~ 3 , as derived from X-ray observations. 
Such an estimate assumes that the hot plasma occupies 
the whole volume of the Coma halo. Assuming that the 
volume where particles are accelerated is V acc = 7.7 x 10 74 
cm 3 , we estimate an optical depth of the acceleration re- 
gion (i.e., a part of the total optical depth where the elec- 
tron spectrum is distorted by the acceleration) which is 
r acc « 9 x 10- 4 . 

The difference between the pure thermal and the total 
(thermal plus quasi-thermal) Sunyaev-Zeldovich effects is 
shown in Fig. 11. This figure shows that the presence of 
sub-relativistic electrons in Coma produces a further tem- 
perature decrement at v ^ 230 GHz and a further tem- 
perature increase at higher frequencies. At the same time, 
the additional pressure (or energy density) contributed by 
the sub-relativistic electrons produce an increase in the 
frequency of the zero point (i.e., the frequency at which 
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Fig. 11. Difference between the pure thermal and the total 
(thermal plus quasi-thermal) Sunyaev-Zeldovich effect. 

AT = 0) of the SZ effect, a value which is uniquely de- 
termined by the overall pressure P to t — Pth + Psubrei , of 
the electron population (see Colafrancesco et al. 2003). 
Also the overall amplitude of the temperature decrement 
due to the subrelativistic electrons is proportional to their 
pressure since AT oc J dlP su brel- 

The change in the Sunyaev-Zel'dovich effect due to the 
accelerated electrons produces a temperature decrement 
AT - 50 - 40 /iK in the frequency range 50 - 100 GHz. 
However, the uncertainty of the available SZ data for 
Coma does not allow yet to set definite constraints on 
the model presented in this paper. Nonetheless, with the 
new generation of telescopes for observing the Sunyaev- 
Zel'dovich effect, such a signal should be detectable. The 
major difficulty in seeing this departure from the ther- 
mal Sunyaev-Zel'dovich effect is the presence of back- 
ground fluctuations in the microwave background radia- 
tion. The expected background anisotropy on degree an- 
gular scales of interest for the Coma cluster is ~ 40 fiK 
in this frequency range, and has a flat spectrum which 
is not readily distinguished from the gently-curving spec- 
trum seen in Figs. 10-11. However, the use of multi fre- 
quency observations increases the possibility of detection 
of the SZ effect associated to sub-relativistic particles since 
the relative amplitudes of the thermal and subrelativis- 
tic SZ effects change with the frequency over the whole 
range (~ 30 — 300 GHz) accessible to SZ experiments. 
Some possibility of detecting the distorted spectrum of the 
Sunyaev-Zel'dovich effect in the presence of background 
fluctuations also exists through a detailed spectral mea- 
surement of the cluster over a smaller region to reduce the 
background fluctuation signal (although this won't reduce 
the amplitude of the kinematic Sunyaev-Zel'dovich effect, 
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which has the same spectrum as the background fluctu- 
ations). Measurements of this type are challenging, but 
should become possible as the sensitivity of bolometer ar- 
rays in the mm and sub-mm bands increases, provided 
that acc urate cross-calibration in the d ifferent bands is 
achieved I Birkinshaw fc Lancaster! 120051) . 

The possibility to measure the SZ effect from the popu- 
lation of sub-relativistic or relativistic particles could pro- 
vide a way to estimate the overall pressure of such particle 
population and, in turn, to constrain the energy spectrum 
taken up by such particles in the cluster atmosphere. Even 
if this demands a large experimental effort, the result will 
undoubtedly shed additional, and maybe crucial, light on 
the nature of the acceleration mechanism. 



8. Conclusions 

We have analyzed in this paper the process of particle ac- 
celeration from a background plasma acting through mag- 
netic fluctuations generated by intraclustcr turbulence 
with a specific application to the problem of the origin 
of hard X-ray emission from the Coma cluster. From the 
equations describing the influence of Coulomb collisions, 
we derived the entire electron distribution function from 
the thermal to the high-energy non-thermal regime. Our 
analysis allowed us to estimate the energy supply neces- 
sary for bremsstrahlung to be responsible for the hard 
X-ray emission. For nonthermal electrons we confirm that 
the bremsstrahlung efficiency is low, which makes it al- 
most impossible to regard this electronic component as 
the source of the hard X-ray excess in Coma. This re- 
sult is in complete agreement with earlier conclusion (see 
Pctrosian 2001) on the inapplicability of the nonther- 
mal bremsstrahlung interpretation. However, we have also 
shown here that the bremsstrahlung efficiency increases 
significantly if the emitting electrons belong to the ex- 
tended transfer regime between the thermal and nonther- 
mal parts of the electron distribution function. This quasi- 
thermal regime is formed naturally when emitting parti- 
cles are accelerated from the background plasma. In the 
specific case of Coma, we found that the total energy loss 
rate of the quasi-thermal electrons that emit the HXR ra- 
diation in the 20 — 80 keV range is almost two orders of 
magnitude lower than for nonthermal particles. This re- 
sult may solve both the problem of the origin of the HXR 
emission in Coma and of the excessive heating of the clus- 
ter gas in the bremsstrahlung interpretation of the HXR 
excess. 

We have further shown that the distribution function 
of quasi-thermal electrons that we derived implies signif- 
icant distortion of the thermal Sunyaev-Zeldovich effect 
from the Coma cluster. Although this additional signal 
is at the level of ~ 10% of the amplitude of the ther- 
mal Sunyaev-Zel'dovich effect and therefore its observa- 
tion will be challenging from the experimental side, its 
definite detection will nonetheless be able to provide a 
stringent test of our theoretical model. 
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Appendix A: kinetic equation 

In the nonrelativistic energy range when E 3> kT (p > 1), 
the kinetic equation for p < p and p > p can be significantly 
simplified and written in the form respectively 



df 
dt 



l_d_ 

p 2 dp 



1 _ 9+2 \ df 

- + ctp q ) + f 

p J dp 



0, 



df 15/15/ 



dt p 2 dp \p dp + ^ 



(A.l) 



(A.2) 



with the boundary conditions corresponding to continuity of 
the function 



f(P,t) = f(p,t), 



(A.3) 



where / is the solution of Eq. 1A.1I for the unknown flux of 
particles, —5*0, running away into the acceleration region 



f{P,t) = \j -n(t) exp 



B(v) 
A(v) 



dv 



(A.4) 



-S 



dv 



exp 



B(t) 
A(t) 



dt 



The flux of particles running into the acceleration region 
derived from Eq. HA .41 is 



S(p) = -So 



v 2 exp (— w 2 /2) dv 



(A.5) 



and this changes from S = for p — to S = — So for p m 
p~M = a~ 1 ^ q+3 ' . Above the momentum p~M, the function / 
cannot be described by a Maxwellian distribution because the 
equilibrium conditions are violated by the run-away flux S, but 
the spectrum is still formed by Coulomb collisions. For energies 



above p~i 



-1/(9+1) 



, the spectrum is nonthermal because 



Coulomb collisions are unimportant particle interactions while 

plasma waves dominate. 

It was shown by IPogieJ l)2000h and iLiane et all <|2002h 

that two excesses above the thermal Maxwellian spectrum are 
formed in the range p > p~M ■ When p~M < p < Pinj , the excess is 
formed by Coulomb collisions (the collisional regime of quasi- 
thermal particles) , and one can imagine the spectrum there as 
a distorted Maxwellian function. For p > p in j the spectrum 
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is formed by particle interactions with plasma waves (the col- 
lisionless regime of nonthermal particles). For these energies 
the particle spectrum can be described as power-law over an 
extended energy range. We also impose a natural boundary 
condition for p — oo, 



/(oo)=0. 



(A.6) 



In order to derive the unknown constant So we should 
match solutions of Eqs. © and l|A,2(l for p — p. 
It is convenient to reduce Eq. )A,2t to the form 



by introducing the variable 



,1/4' 



5^2 
and the function z 




V2 9 / 4 



4/5 



(A.7) 



(A.8) 



(A.9) 



(A.10) 



(A.11) 



For large £ (or equivalent p or £) values we can neglect the 
right-hand side in Eq. 1A.7L and the equation takes the well- 
known one-dimensional diffusion equation form 



a"z 



dz 
dt 



0, (A.12) 
with boundary condition JA.3I . As an initial condition we put 



5(0 =0 



at 



t = 0, 



(A.13) 



i.e., we assume that there were no particles in the momentum 
range p > p at t = 0. 

The solution of Eq. lEOl can be presented with the well- 
known Green function G(x, xq \ t, r) for the diffusion equation 



G(x,x I t,r) 
1 1 



(A. 14) 



20F (t 



\l/2 



exp 



(x — x ) 2 



4(t-r) 



Then the function / for Eq. HA.12I with the boundary condi- 
tions JA.3t has the form 



/(P,*) = /(P) 



1 - $ 



p exp(p 2 ) 
p exp(p 2 ) 

(C - 

2^/t 



1/4 



(A.15) 



where Q(x) is the error function. We see that for t — > 00 the so- 
lution (IA. 1511 tends to the stationary distribution. In this limit 
df /dt <C 1 and our neglect of the right-hand side of Eq. I A. 71 
is not valid. Therefore in the stationary distribution limit, the 



distribution function / tends to the equilibrium Maxwellian 
function for p > p when t — > 00, 



f(p, i = 00) = C(t) exp 



P_ 
' 2 



(A.16) 



The constant C and the unknown flux So can be defined 
from the boundary condition 1A.3I and the normalization con- 
dition 



f(P)p dp — Sot. 



(A.17) 



However, as we see from Eq. 1A.15L in the stationary case we 
have an exponential cut-off of the distribution function / for 
p = p, which gives an approximate equation for So of the form 



So = \l -n(t) 



rp dv 
x I / , ; , exp 



(A.18) 



B(t) 
A(t) 



dt 



where n(t) is the density of background plasma, which de- 
creases slowly with time. 
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